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Introduction 


As  a  part  of  the  problem  of  developing  numerical  wave 
forecasting  procedures  for  the  North  Atlantic  Ocean,  selected  se¬ 
quences  of  the  weather  maps  for  the  North  Atlantic  for  which  wave 
data  were  known  to  be  available  were  studied  in  detail  for  the  five 
year  period  beginning  in  April  1955  and  ending  in  March  I960. 

Certain  dates  and  times  of  observations  were  selected  for  a  variety 
of  reasons  for  study.  For  these  dates  and  times,  the  National  Insti¬ 
tute  of  Oceanography  provided  copies  of  the  wave  records  that  were 
obtained  by  the  OWS  Weather  Explorer  and  by  the  OWS  Weather  Reporter 

In  total,  about  800  wave  records  were  provided,  and  a 
complete  spectral  analysis  is  planned  for  about  400  of  these  records. 

This  report  is  the  first  of  a  series  of  reports  to  present 
in  tabular  and  graphical  form  the  results  of  these  analyses.  The  total 
number  of  spectra  given  is  114. 

Analysis  procedures 

The  original  wave  records  varied  in  length,  but  almost  all 
of  this  first  set  were  15  minutes  long.  The  crest  to  trough  heights 
of  the  highest  waves  in  a  particular  record  (uncorrected  for  cali¬ 
bration  effects)  varied  from  a  few  feet  to  more  than  60  feet  in  the 
complete  set  of  records.  Bounds  were  set  on  each  record  just  above 
the  highest  wave  crest  and  just  below  the  lowest  wave  trough,  and 
the  records  were  read  to  an  accuracy  of  one  part  in  a  thousand 
(nominally)  over  this  range  at  an  interval  of  1.5  seconds  throughout 
the  record.  Thus  a  15  minute  record  was  reduced  to  a  time  series 


of  600  points.  Where  gaps  or  irregularities  occurred,  the  records 
were  smoothed  by  hand  as  accurately  as  possible. 

The  time  series  of  600  points  was  then  analyzed  on  the 
CDC  1 6 0 4  so  as  to  estimate  the  energy  spectrum  of  the  waves  at  60 
points  over  the  frequency  range  from  zero  to  0.333  cycles  per  second 
by  means  of  the  procedures  given  by  Tukey  (1949)  as  explained  in 
detail  by  Blackman  and  Tukey  (1  958).  The  smoothing  operation  that 
was  used  to  go  from  L  to  U  in  the  equations  of  Blackman  and  Tukey 
was 

'  h  h- 1  h  h+1 

with  suitable  corrections  at  the  ends  of  the  range. 

The  spectral  estimates  so  obtained  still  had  to  be  cor¬ 
rected  for  the  response  of  the  shipborne  wave  recorder  (Tucker, 

1956)  and  for  the  introduction  of  noise  in  both  the  original  record 
and  in  the  digitization  procedure.  The  calibration  of  the  shipborne 
recorder  depends  on  the  ship,  and  the  calibration  curves  were  pro¬ 
vided  by  Mr.  D.  E.  Cartwright  for  this  purpose.  The  calibration  curve 
for  the  Weather  Explorer  is  given  by  Table  1.  The  calibration  curve 
for  the  Weather  Reporter  is  given  by  Table  2. 

As  in  another  investigation  (Bretschneider ,  Crutcher,  et 
al  (in  press)),  it  was  found  that  the  application  of  the  above  cali¬ 
bration  curves  to  the  spectra  that  were  originally  computed  resulted 
first  in  a  decrease  and  then  a  rapid  increase  in  the  spectra  at  high 
frequencies  due  to  the  presence  of  noise  and  other  irregularities 
(possibly  from  nonlinear  effects  in  the  original  wave  records)  at 
the  high  frequency  end  of  the  spectrum.  To  eliminate  this  effect, 
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the  last  part  of  the  spectrum  was  smoothed  by  a  three  point  running 
weighted  mean  (0.25,  0.50,  and  0.25)  and  then  the  last  ten  values 
were  averaged.*  This  average  was  treated  as  white  noise  and  sub¬ 
tracted  from  all  spectral  estimates.  When  the  reduced  values  were 
multiplied  by  the  appropriate  calibration  curves,  the  usual  result 
was  a  fairly  smooth  spectrum  that  decreased  regularly  toward  zero 
values  at  high  frequency.  By  such  a  procedure  some  of  the  spectral 
values  at  high  frequency  will  be  negative.  These  values  were  auto¬ 
matically  set  equal  to  zero  in  the  rest  of  the  computations. 

Inadvertently,  Table  1  was  applied  where  Table  2  ought 
to  have  been  applied  to  some  of  the  spectra.  The  result  was  ex¬ 
ceptionally  high  values  at  high  frequencies.  Thus  the  calibrations 
given  in  these  two  tables  do  distinguish  between  the  characteristics 
of  the  two  ships  as,  when  done  correctly,  reasonable  results  were 
obtained. 

Even  with  these  corrections,  there  were  a  few  spectra  that 
still  became  exceptionally  large  for  frequencies  greater  than  about 
0.25  cycles  per  second.  This  behavior  was  apparently  caused  by  the 
original  quality  of  the  record  and  not  by  the  digitization  procedure. 
These  spectra  were  further  modified  by  arbitrarily  setting  the  cali¬ 
bration  curve  equal  to  one  above  a  certain  frequency  that  was 
selected  by  inspection  of  each  spectrum. 

The  result  of  such  a  sequence  of  computations  should  yield 
fairly  reliable  spectral  estimates  for  frequencies  ranging  from  zero 
to  0.25  cycles  per  second,  but  the  values  at  high  frequencies  should 
not  be  used  to  decide  on  any  features  of  the  high  frequency  end. 


*See  p.  6  for  details. 
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Sample  parameter  estimates 

The  spectral  estimates  that  resulted  from  this  sequence  of 
operations  were  then  processed  further  to  obtain  some  additional 
useful  information.  Let  for  h  =  0,  1,  2,  •  •  •  ,  60,  represent  the 

spectral  estimates  (after  subtraction  of  the  noise  and  multiplication  by 
the  calibration  for  the  shipborne  recorder)  in  terms  of  the  resolution 
of  the  variance  of  the  wave  record  into  frequency  intervals.  The 
following  quantities  were  then  also  computed  and  tabulated  with  each 
spectrum. 

(2)  CORR  VAR  =  corrected  variance  = 

(3)  STG  HGT  =  hL  =  2.83  (2ZU,  *)1/2 

a  n 

~  1/2 

(4)  AVER  T  =  T  =  [SU^'/Sf^U^] 

(5)  TOTAL  DF  =  Total  degrees  of  freedom  =  10  [2  Uh*]2/  [S  ] 

(for  600  points,  60  lags;  i.  e.  ,  20  degrees  of  freedom 
per  spectral  estimate) 

The  confidence  intervals  on  the  corrected  variance  and  on  the 
significant  height  are  given  by 

Upper  95%  on  CORR  VAR  =  )  CORR  VAR 

(6) 

Lower  5%  on  CORR  VAR  =  ( 1  0  7  )  CORR  VAR 


and  by 


(?) 


rr  ncrf  74  1a+1/2\/TDF  75 

Upper  95%  on  =  10  Hi 

3  3 

Lower  5%  on  Hi  =  10'1/2n^°"^Hi 

3  3 

in  terms  of  the  total  degrees  of  freedom  (TDF)  to  a  high  degree  of 
accuracy  since  the  total  degrees  of  freedom  are  large. 

The  corrected  variance,  the  significant  height,  and  the 
total  degrees  of  freedom  are  relatively  insensitive  to  changes  in  the 
noise  level  and  in  the  high  frequency  behavior  of  the  spectrum. 
However,  the  average  period  can  properly  be  viewed  with  caution. 

The  winds  near  the  ship  at  the  time  of  observation  are 
also  given  to  the  nearest  five  knots  as  read  directly  from  weather 
maps.  These  values  are  subject  to  later  correction  in  terms  of  the 
logs  of  the  weather  ships. 

Explanation  of  tables  and  graphs 

The  body  of  this  report  consists  of  supplementary  tables, 
of  tables  that  give  the  appropriate  results  for  each  of  the  original 
wave  records,  and  of  graphs  of  each  of  the  estimated  spectra  along 
with  the  confident  intervals  on  the  spectra. 

The  supplementary  tables  consist  of  Tables  1  through  4. 
Tables  1  and  2  have  been  described  above. 

Table  3  gives  either  the  on  station  position  of  the  ship, 

A,  I,  J,  or  K,  or  the  latitude  and  longitude  of  the  ship  if  it  is  going 
on  or  off  station.  The  speed  and  direction  of  the  ship  is  given. 


Position  A  corresponds  to  62  “N,  33  °W. 

Position  I  corresponds  to  59  #N,  19  "W, 

Position  J  corresponds  to  52.5°N,  20°W. 

Position  K  corresponds  to  45  °N,  16°W. 

If  the  record  was  not  15  minutes  long,  less  than  600  points  were 
read.  For  these  records,  Table  4  gives  the  actual  number  of  points 
used  and  the  corrected  total  degrees  of  freedom.  A  correction  to  the 
upper  and  lower  confidence  limits,  which  would  be  quite  small,  would 
also  be  needed  to  be  exact. 

Spectral  tabulations 

A  tabulated  spectrum  can  be  interpreted  as  follows: 

la)  Supplementary  data  for  each  spectrum  consist  of  the  date, 

hour,  wind  speed,  total  degrees  of  freedom,  average 
period,  significant  height,  corrected  variance,  noise 
level,  and  record  number.  Some  tables  give  the  confi¬ 
dence  limits  for  the  height  according  to  equation  (7). 

lb)  In  the  first  column,  the  spectral  lag  numbers  (H)  are  given. 

2)  In  the  second  column  (FRE)  the  frequency  according  to  the 

equation  f  =  H/180  (sec  ')  is  given. 

2 

3)  In  the  third  column  (UNIT  =  FT  ),  the  spectrum  as  com- 

2 

puted  from  the  original  data  is  given  in  units  of  (ft)  . 

4)  In  the  fourth  column  (FILTERED)  a  smoothing  operator 

for  H  >  40  is  applied.  It  is  actually 

Fh  =  0.25Uh_j  +  0.50Uh  +  0.25Fh+1 
(where  F  =  Filtered,  and  U  =  Unit) 

5)  In  the  fifth  column  (LESS  NOISE),  the  noise  level  shown  at 

the  top  is  subtracted  from  each  estimate. 
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6)  In  the  sixth  column  (CORR  FT  2),  the  LESS  NOISE 
column  is  multiplied  by  the  calibration  curve  for  the 
shipborne  record  according  to  either  Table  1  or 

Table  2.  If  this  column  agrees  with  the  previous  column, 
at  high  frequencies,  the  calibration  curve  has  been  arbi¬ 
trarily  set  equal  to  one  to  avoid  extreme  values  at  high 
frequency. 

7)  In  the  last  two  columns,  the  upper  and  lower  95%  and  5% 
confidence  bounds  are  shown. 

The  graphs  of  the  spectra 

The  graph  that  accompanies  the  spectral  tabulation  shows  the 

spectrum  and  the  95%  and  5%  confidence  bounds.  The  scale  is  chosen 

so  that  the  highest  95%  confidence  value  is  at  the  top  of  the  graph  and 

the  vertical  axis  of  the  coordinate  system  shows  the  spectral  values 

2 

for  that  spectrum  in  units  of  (feet)  .  The  scales  change  with  each 
spectrum  and  comparisons  between  spectra  by  means  of  the  graphs 
should  be  made  cautiously. 
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Table  1.  Calibration  factors  for  the  Weather  Explorer  . 


1.0000 

1.0C00 

1.0000 

1. 

OCOO 

1.0000 

1.0000 

1.6157 

1.3740 

1 .2452 

1. 

1746 

1.1399 

1.1291 

1.1343 

1.1547 

1.1870 

1. 

2304 

1.2845 

1.3504 

1.4277 

1.5193 

1.6241 

1. 

7444 

1.8828 

2.0415 

2.2243 

2.4349 

2.6765 

2. 

9523 

3.2725 

3.6414 

4.0714 

4.5654 

5.1490 

5. 

8190 

6.0136 

7.5383 

8.633b 

9.9169 

1 1.4459 

13. 

2691 

15.4245 

18.U095 

21. 1086 

24.8366 

29.3522 

34. 

8079 

4 1.4485 

49.5464 

59.4548 

71.5502 

86.5947 

105. 

1503 

128.1186 

156.7723 

192.5202 

237.3987 

293. 8682 

365  . 

1736 

455.5306 

570.2699 

710.87C5 

Table  2.  Calibration  factors  for  the  Weather  Reporter. 


1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.5755 

1.3277 

1 . 19C8 

l .  1 C99 

1.0630 

1.0375 

1.0257 

1.0260 

1  .0350 

1.0514 

1.0633 

1 .1034 

1.1384 

1 .1805 

1.2280 

1.2817 

1.3424 

1.4105 

1.4871 

1.5731 

1  .6684 

1.7736 

1.8918 

2.0229 

2.1704 

2.3321 

2.5169 

2.7161 

2.9479 

3.2018 

3.4899 

3.8C88 

4.1715 

4.5826 

5.0408 

5.5616 

6.1512 

6.8201 

7.5845 

8.4517 

9.4439 

10.5785 

11.8784 

13.3689 

15.0856 

17.0596 

19.3530 

22.0055 

25.0761 

28.6529 

32.8206 

37.6868 

43.3807 

50.0432 

57.8872 

Table  3.  Position  and  speed  of  ship  for  each  record 


Record  No. 

DL  1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 


Position 

Heading 

Speed  (kts) 

I 

stopped 

I 

stopped 

I 

stopped 

I 

stopped 

I 

stopped 

I 

stopped 

I 

stopped 

I 

stopped 

I 

345  0 

1 

I 

360° 

2 

I 

360° 

1 

I 

stopped 

I 

stopped 

I 

stopped 

J 

280  ° 

1 

J 

280  ° 

1 

J 

stopped 

J 

290° 

1 

J 

290 0 

1/2 

J 

stopped 

J 

stopped 

J 

280° 

1 

J 

280° 

1 

J 

hove  to 

J 

stopped 

J 

260° 

1/2 

J 

280° 

0 

J 

280° 

1/2 

J 

stopped 

J 

stopped 

Table  3.  (cont.  ) 


Record  No.  Position  Heading  Speed  (kts) 


31 

I 

stopped 

32 

I 

stopped 

33 

I 

190° 

1 

34 

I 

vO 

o 

o 

2 

35 

I 

250° 

1 

36 

I 

230° 

1 

37 

I 

230° 

1 

38 

I 

240° 

1/2 

39 

I 

240° 

1/2 

40 

I 

240° 

1 

41 

I 

250° 

1 

42 

I 

stopped 

43 

A 

070° 

2 

44 

A 

stopped 

45 

A 

stopped 

46 

I 

250° 

1 

47 

I 

250° 

1} 

48 

I 

250  ° 

2 

49 

I 

250° 

2 

50 

I 

250  ° 

1 

51 

I 

260° 

1 

52 

I 

265  ° 

4 

53 

I 

260° 

2 

54 

I 

265  ° 

1 

55 

I 

280° 

2 

56 

I 

stopped 

57 

I 

305  ° 

4 

58 

I 

stopped 

59 

I 

stopped 

60 

I 

stopped 

Table  3.  (cont.  ) 


Record  No. 

Position 

Heading 

Speed  (kts) 

rH  1 

A 

040° 

1 

2 

A 

090 6 

2 

3 

A 

050° 

2 

4 

A 

040° 

1 

5 

A 

040° 

1 

6 

A 

stopped 

7 

A 

stopped 

8 

K 

stopped 

9 

K 

stopped 

10 

K 

stopped 

11 

K 

290° 

1 

12 

K 

280° 

1/2 

13 

K 

275  ° 

1/2 

14 

K 

275° 

1 

15 

K 

275  ° 

4 

16 

K 

270° 

1 

17 

K 

275  ° 

1 

18 

K 

275  ° 

1 

19 

K 

stopped 

20 

A 

085  ° 

2 

21 

A 

140° 

2 

22 

A 

160° 

23 

A 

240  ° 

1 

24 

A 

200° 

2 

25 

A 

200° 

2 

26 

A 

220° 

2 

27 

A 

250° 

1 

28 

A 

240° 

2 

29 

A 

235  ° 

30 

A 

240° 

2 
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Table  3.  (cont.  ) 


Record  No. 

Position 

Heading 

Speed  (kts) 

JH  31 

A 

230° 

32 

A 

190° 

2 

33 

A 

190° 

2 

34 

A 

o 

O 

oo 

1— t 

2 

35 

A 

s  topped 

36 

A 

230° 

37 

A 

240° 

2 

38 

A 

230° 

la1 

39 

A 

230° 

2 

40 

A 

220° 

1 

41 

J 

080° 

1 

42 

J 

090° 

1 

43 

J 

070° 

1 

44 

52  “42  'N,  1 9  ° W 

065° 

3 

45 

53  °N  ,  1 8 0 W 

060° 

7 

46 

52  °54'N,  16  °54' W 

060° 

7 

47 

53  ”01  'N,  16  °36'W 

050° 

4 

48 

53  o06'N,  16  °20'W 

350° 

4 

49 

53  °08'N,  16’15'W 

050° 

5 

50 

53  0  1  8'N,  16  °1 6' W 

340° 

1 

51 

53  °  1  8'N,  1 6 °  1  6 1 W 

hove  to 

52 

53  ° 1 7 'N, 1 5  °56' W 

120° 

1  1 

53 

52  0 48'N,  14°26'W 

120° 

10 

54 

52  °30'N,  13°35'W 

120° 

Hi 

Table  4.  Data  on  short  records  for  which  less  than 
600  points  were  available. 


Record  No. 

no.  of 
points 

Original 

TDF 

Corrected 
_ TDF 

DL  50 

592 

144 

142 

JH  4 

561 

150 

140 

JH  16 

591 

134 

132 

JH  17 

592 

101 

100 

JH  18 

590 

128 

126 

JH  19 

581 

133 

129 

JH  24 

5  86 

151 

147 

JH  36 

585 
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199 
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